Pushing off the walls: a mechanism of cell motility in confinement 
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We propose a novel, mechanism of cell motility, which relies on the coupling of actin polymer- 
ization at the cell membrane to geometric confinement. We consider a polymerizing viscoelastic 
cytoskeletal gel confined in a narrow channel, and show analytically that spontaneous motion oc- 
curs. Interestingly, this does not require specific adhesion with the channel walls, and yields velocities 
potentially larger than the polymerization velocity. The contractile activity of myosin motors is not 
necessary to trigger motility in this mechanism, but is shown quantitatively to increase the veloc- 
ity. Our model qualitatively accounts for recent experiments which show that cells without specific 
adhesion proteins are motile only in confined environments while they are unable to move on a flat 
surface, and could help in understanding the mechanisms of cell migration in more complex confined 
geometries such as living tissues. 

PACS numbers: 87.10.-e, 87.17.Jj, 83.80.Lz 



Besides its obvious importance for biology, cell motility 
has motivated numerous studies in the physics commu- 
nity. Identifying simple mechanisms of self-motion of soft 
condensed matter is an important challenge for physics, 
cell biology and biomimetic material technology. Sus- 
tained motion at low reynolds number necessitates a con- 
stant energy input, and therefore an active system, that 
is a system driven out of equilibrium by an internal or an 
external energy source. The cell cytoskeleton is a strik- 
ing example of such an active system. It is a network of 
long semi-flexible filaments made up of protein subunits, 
interacting with other proteins such as motor proteins 
which use the chemical energy of ATP hydrolysis to ex- 
ert active stresses that deform the network [1]. Other 
examples of active systems range from animal flocks to 
bacterial colonies [3,11] and vibrated granular media 

Modelling cell motility has inspired both experimen- 
talists and theoreticians, who have now distinguished 
two main different mechanisms: polymerization (tread- 
milling) and contractility. Both polymerization induced 
motion Q and contractility induced spontaneous flows 
[l[ have now been observed in- vitro, and studied theoret- 
ically @, @, and numerically [l(| • In all models the key 
ingredients of motion are an energy input to compensate 
dissipation and sufficient adhesion or friction with a sub- 
strate to acquire momentum. The usual picture of cell lo- 
comotion is then as follows: the cell lamcllipodium builds 
strong adhesion points with the substrate and pushes for- 
ward its membrane by polymerizing actin. At the back, 
the cell body contracts and breaks the adhesion points. 
In particular, the overall cell velocity is then limited by 
the actin polymerization rate (which however varies sub- 
stantially between cell types) , in agreement with avail- 
able experimental dat a fill . Ha . HH . 

In a recent paper 14 1 (see also [l5[ for another cell 



type), it has been observed both in- vivo and in- vitro 
that mutant dendritic cells (DC) that are unable to pro- 
duce active integrin complexes (adhesion proteins) dis- 
play sustained motility in confined environments (tissues 
or synthetic polymeric gels), whereas they fail to move 
on flat two-dimensional substrates due to their reduced 
adhesion ability. These observations suggest the exis- 
tence of an alternative mechanism of motility to the ad- 
hesion dependent picture outlined above. Here we pro- 
pose a new, simple mechanism of motility which accounts 
for these observations. This mechanism is mainly pow- 
ered by actin polymerization at the cell membrane, and 
strongly relies on geometric confinement. Interestingly, it 
does not necessitate strong specific adhesion, and yields 
velocities potentially larger than the polymerization ve- 
locity. This confinement induced motility mechanism is 
backed by in- vitro experiments of DC motility in micro- 
fabricated channels [16| (see Fig. [I]). 

We first introduce our model in its minimal form of a 
polymerizing visco-elastic gel confined in a channel. We 
then refine this model to mimic motile cells in confine- 
ment. Finally, we show quantitatively that the contrac- 
tile activity of myosin increases the velocity of motion. 

The model, which relies on the hydrodynamic theory of 
active gels 0, is as follows. We consider an incompress- 
ible viscoelastic film confined in a bidimensional channel 
of width b. Note that this bidimensional geometry mim- 
ics the experimental conditions of channels of rectangular 
section (Fig. [I}, and that the case of a cylindrical con- 
fining channel can be treated with minor modifications. 
The axes are defined with x along the channel and z 
across it. The confining walls are placed at z — and 
z = b. We denote the components of the velocity of the 
gel by «i, and the strain rate by ity = (diVj + djVi)/2. 
We assume that the gel is described by a linear Maxwell 
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Defining the average velocity along the channel v(x) — 
(1/6) L v x (x, z)dz, we obtain the following Darcy's law: 
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FIG. 1: a) RICM image of a dendritic cell moving (to the 
right) in a channel of 4/jm width. The dark zone at the back of 
the cell (left) indicates a large contact of the membrane with 
the channel wall (independent of the nucleus, dotted line), 
and therefore a high normal constraint, compared to the front 
(right). The typical observed velocity reaches 12 — 15/^m/min 
in channels and 4 — 6/^m/min on a flat surface, b) Channel 
geometry and model. The arrows show the flow direction. 



model of visco-elasticity. The constitutive equation re- 
lating the strain rate to the deviatory stress tensor cr^ is 
then written 2rjUij = (l + t^) ov, , where r\ is the shear 
viscosity, r is a typical relaxation time and D/Di de- 
notes the convective derivative. For dilute polymer gels, 
r is very small and the gel behaves as a viscous fluid. 
At higher concentrations, or for more cross- linked gels, 
r becomes very large and the gel behaves as an elastic 
medium. In what follows we assume that the gel is in 
either of the two regimes. For the sake of simplicity we 
assume here that the gel is incompressible, which means 
that r depends on pressure P only. We define a criti- 
cal pressure such that t(P < P*) = (viscous regime) 
and t(P > P*) — oo (elastic regime). We next suppose 
that the gel is polymerized at the gel/substrate interface 
with speed v p = v z (x,z = 0) = — v z (x,z = 6) in the 
viscous regime, as depicted in Fig. [TJ This is justified by 
the common observation of actin polymerization activa- 
tors such as WASP proteins preferentially located along 
the cell membrane [ill E3- I n the case of DCs, actin 
filaments can anchor perpendicularly to the cell mem- 
brane, forming structures called podosomes where poly- 
merization takes place, therefore inducing an inward flow 
of actin as it is assumed in our model [181 ]. Finally we 
assume viscous friction at the channel walls z = and 
z = b, and write o~ xz = £v x where £ characterizes the 
friction. 

We now derive the dynamical equations of the sys- 
tem in the lubrication approximation where L 
is the typical length of the system). In this limit the 
Reynolds number is small and the velocity field Vi(x, z) 
can be obtained from the force balance d x P = T]d^v x 
and the condition of incompressibility u xx + u zz = 0. 



Including the depolymerization of the gel kd, mass con- 
servation of the gel reads ^ = 2v p /b — kd, which gives 
in turn: 
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where v p and the nondimensional friction £ = £6/677 can 
be a priori functions of P and x. 

Two boundary conditions are needed to determine the 
pressure profile P(x). We neglect the friction with the 
surrounding fluid in the channel and set the pressure at 
the leading edge, which is assumed to coincide with the 
point x = L, as P(L) = 0, which gives the first bound- 
ary condition. Note that if the pressure at the leading 
edge is finite due to an external force, our results ap- 
ply with an unimportant shift in the pressure field. We 
look for stationary states with broken symmetry and pos- 
itive velocity and therefore the pressure is a decreasing 
function of x. We then argue that if the system is large 
enough, there exists a travelling front of gel of length L 
in the fluid phase, travelling at velocity V. The back 
boundary of this front fluid part coincides with the point 
x = where the pressure reaches the threshold P* , be- 
hind which is a growing elastic part. Such a denser elastic 
region at the back of DCs, called the uropod, is indeed 
well reported and is characterized by a higher con- 
centration of cross-linkers. As the velocity of the elastic 
part should be zero one has v(0) = = ^\ x =o, giving 
the second boundary condition which allows the explicit 
calculation of the pressure field. The self-consistent con- 
dition P(0) = P* gives in turn an equation enabling the 
calculation of the length L of the fluid front. We then 
write that the velocity V of the front is given by the 
calculated velocity of the flow plus the polymerization 
velocity at the leading edge v{L) + v p (L). We stress that 
the flow velocity is forward, i.e. in the same direction 
as the moving leading edge. Note that the length L of 
the fluid front is constant, which necessitates that the 
elastic/fluid boundary moves at the same velocity V. 

Qualitatively, the value of the length L is dictated by 
the steepness of the pressure gradient, and therefore by 
the friction £. If £ is small, then only very long fluid 
fronts can move. We show now quantitatively that the 
coupling of £ with the pressure field actually enables short 
fluid fronts to move even with a low bare friction. The 
key ingredients are as follows. Following [2(| we argue 
that the friction coefficient £ depends on the normal con- 
straint in the case of a polymeric gel. Indeed qualitatively 
a high normal constraint increases the attachment rate 
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of polymers onto the channel walls by lowering the en- 
tropic barrier, and decreases the detatchment rate. It is 
shown in [20| that in the regime of moderate tangential 
speed, one has £ = £ t0 e^ tyP ~' Jnn \ where in our geometry 
the normal stress o~ nn = a zz for both walls at z = and 
z = b. Next, following standard ratchet models 21 1, the 



polymerization speed at the cell membrane is assumed to 
depend on the normal constraint in the gel according to 
v p = v < p ] e~ a ( p ~ cr " s> . These assumptions make Eq. ([2]) an 
autonomous equation for P, which is completed by the 
two boundary conditions discussed above. 

To our knowledge, such an equation cannot be solved 
analytically in the general case. In the regime of small a 
(defined by aP* <C 1), v p can be taken as constant and an 
analytical approximation scheme can be proposed, which 
enables a discussion of the motility mechanism. We first 
neglect the elongational shear stress and write £ « £,oe ■ 
This assumption underestimates the friction, and there- 
fore the pressure field. Eq. [2] can then be integrated and 
yields an implicit equation for the pressure field to lowest 
order P°: 
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This first expression P° gives a lower bound of the 
pressure field, and provides a satisfactory approximate 
(Fig. [2|). To go further, we use P°(x) to determine the 
lowest order velocity profile v°(z) and calculate a° z — 
2ijd z v®. Eq. ^ is then resolved using this calculated nu- 
merical value of a® z , yielding the next order P 1 (x) which 
in turn can be used for further iterations. Using realistic 
values for the parameters corresponding to the actin cy- 
toskeleton, used in [2(| and [2§| , the procedure converges 
rapidly. Note that in the general case of any a, the first 
iteration giving P°(x) has to be performed numerically, 
and then the same procedure applies. 

This mechanism therefore produces a forward flow 
which relies on a pressure build up to P* in the gel, here 
induced by confinement. It requires a minimal system 
size L given by taking x — in Eq. [3J Interestingly, this 
can be obtained even for a very low bare friction coeffi- 
cient £o: since the exponential dependence of the friction 
on the pressure field permits the effective friction to reach 
large values for finite L, enabling motion. 

On the other hand, Eq. [3] shows that L increases 
faster than linearly with b, indicating that this mecha- 
nism would not be significant in the case of a gel on a 
flat open substrate. In this case, which depicts the lamel- 
lipodium of a cell lying on a flat substrate (see p}), the 
typical confining length b is large (of the order of the 
cell size), and the typical length L necessary to build 
a strong pressure gradient is very large (L > cell size). 
As the pressure gradient in the cell is then much weaker 
than in the confined case, the friction with the substrate 
remains close to its bare value, yielding a much smaller 
momentum transfer with the substrate. Additionally, we 
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FIG. 2: Pressure profile. Dashed red: numerical value P° for 
(3 = 1 and a = O.OlkPa -1 , dotted blue: analytical expression 
Eq. (|3]) for P = 1, a = 0, dot-dashed green: numerical value 
P 1 for /3 = 1, a = 0, solid black: numerical value P° for 
P = 0, a = O.OlkPa -1 . Other parameters (estimates from 
[13, [3]): L = Wfj,m, b = 1/xm, rj = lOkPas, k d = 0.1s" 1 , 



= 0.1/ims , £0 = O.lkPa s /jm 



is taken as the 



speed of DCs on a surface which is expected to be the actin 
polymerisation speed [l3[ . £0 is taken as very small (lowest 
estimate in [2(j which is 100 fold smaller than in keratocytes 
pi) to mimic the low adhesion of integrin knockout DCs. 



then expect that the pressure remains below P*, and 
that no elastic phase is formed, thus preventing forward 
flow. Without confinement, our model therefore suggests 
a retrograde flow, that is in the opposite direction to the 
moving leading e dge, a s previously modeled and observed 
for lamellipodia [7|,[23j on flat substrates. The flow direc- 
tion, and therefore the direction of the pressure gradient 
in the gel, constitutes the main difference between the 
confinement induced mechanism of motility that we re- 
port here and the standard picture of cells lying on flat 
substrates. 

Experimentally the pressure field can be quantified in- 
directly by measuring the effective contact area of the 
cell membrane using Reflection Interference Contrast Mi- 
croscopy (RICM). Fig. JlBl shows clearly that for a 
DC confined in a channel a larger contact area at the back 
of the cell is seen, indicating a backward pressure gradient 
in qualitative agreement with our theoretical prediction. 

We now argue that this mechanism of confinement in- 
duced motility could be used by cells such as DCs to 
move in confined environments like channels. Extra hy- 
potheses have to be added to the above model in order to 
more realistically capture the geometry of a moving cell. 
Instead of an open system, we assume now that the back 
edge is a thin slice where the gel is in its elastic regime, 
which mimics the uropod observed at the back of the cell. 
Additionally, we assume (3P* 3> 1 such that the friction 
of the uropod with the channel walls is very large, enforc- 
ing v(0) = 0. To conserve the total cell mass, we further 
assume that in the uropod the gel depolymerizes at the 
speed of the leading edge V = v(L)+v p (L) (which defines 
the over all speed of the cell). A high depolymerisation 
rate in the uropod can be justified by the high pressure 
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and a depletion of free actin monomers due to the forward 
flow. With these hypotheses, the model presented above 
mimics a cell moving in a channel with velocity V, and 
shows that the confinement induced motility mechanism 
can indeed be used by cells. Interestingly, for the param- 
eter values used in Fig. [5] the front velocity is calculated 
to be ~ 10/im/min which is close to the velocity that 
can be reached by DCs in collagen matrices [IH and in 
channels (up to 12 — 15/im/min) id] which significantly 
is larger than the polymeristation velocity taken as the 
speed on a flat surface, 4 — 6/xm/min 16|, 24 1. 



is well captured by our model of confinement induced 
motility. Importantly, this mechanism is widely indepen- 
dent of adhesion properties with the substrate, since the 
mechanism relies on an enhancement of friction due to 
a pressure build-up, and does not require specific adhe- 
sion proteins. In particular, this result is compatible with 
the experiments of 14j, where it is found that integrin 
knocked out DCs are motile only in confined environ- 
ments. Finally, the effect of myosin induced contractility 
can be taken into account, and yields a further enhance- 
ment of motility, in agreement with experiments (13 , [IB] • 
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FIG. 3: Polarization (black) and flow velocity with (dashed 
green or gray) and without (blue or dark gray) myosin. Polar- 
ization arrows point in the direction of actin polymerisation. 
The active term is taken as £(x)A/i = £A[i(L—x)/L for x > 
such that there are more myosins at the back of the cell. 

Finally, we show that the coupling of the contractile 
effect of myosins with the normalized polarization field 
Pi of actin filaments (parametrized by its angle with the 
x axis) can also be taken into account. The polarization 
9{z) = 9 + 8i where 9q = — f (1 — is the static con- 
figuration which satisfies the normal anchoring boundary 
condition and 6\ is assumed to be linear in v®. Taking 
into account the active stress af!f tlvc = ((x)Afj,piPj from 
0, we obtain a perturbative solution around v® = for 
the polarization and flow fields (shown in Fig. [3]). Here 
it is also useful to consider the lubrication approximation 
b — > 0. To lowest order in this approximation, the po- 
larization is given by its static configuration 9 and one 
obtains a generalized Darcy's law: 
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where C(a;)A/i stands for the active coupling of myosins 
to actin filaments (see Q for review), which is to linear 
order proportional to the myosin concentration. This 
equation shows that the contractile active stress induced 
by myosins (£A/z < 0), increases the velocity of the actin 
flow, as shown in Fig. [3] 

In conclusion, the motility mechanism of DCs in con- 
fined environments is strikingly different from the stan- 
dard picture of cell motility on open flat substrates, and 
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